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NEW LOOK AT THE NAVIER-STOKES EQUATION 


TOMASZ DLOTKO 


Abstract. We propose a new way of looking at the Navier-Stokes equation (N-S) in dimensions two 
and three. In 2-D that problem is critical with respect to the standard L 2 a priori estimates. We 
consider its regular approximations in which the —PA operator is replaced with the fractional power 
(-PA) 1 ^, o: > 0 small, where P is the projector on the space of divergence-free functions. The 3-D N-S 
equation is super-critical with respect to the standard L 2 a priori estimates; the regular approximating 
problem in 3-D should contain fractional power (— PA) S with s > 

Using Dan Henry’s semigroup approach we construct regular solutions to such approximations. The 
solutions are unique, smooth and regularized through the equation in time. Solution to 2-D and 3-D N-S 
equations are obtained next as a limit of the regular solutions of the above approximations. Moreover, 
since the nonlinearity of the N-S equation is of quadratic type , the solutions corresponding to small 
initial data and small / are shown to be global in time and regular. 


Dedicated to my Professors: Andrzej Lasotaf Jack K. Ha+ and Nick D. Alikakos. 

1. Introduction 

Fifteen years ago, in the monograph [3], we were studying a direct generalization of semilinear para¬ 
bolic equations, namely abstract semilinear equations with sectorial positive operator in the main part. 
Following the idea of Dan Henry [20], such equations will be studied using classical techniques of the 
theory of ordinary differential equations, modified however to cover equations with unbounded operator 
in a Banach space. Such approach, located inside the semigroup theory, proves its strength and utility 
in the study of several classical problems; some of them were reported in [2BGS- But such technique 
offers further possible generalizations, first to study the problems, like e.g. Korteweg-de Vries equation 
and its extensions [a 0 na os], where the solutions are obtained as a limit of solutions to parabolic 
regularizations of such equations (the method known as vanishing viscosity technique, originated by E. 
Hopf, O.A. Oleinik, P.D. Lax in 1950th). Another possible application of Henry’s technique is to study 
critical problems (e.g. mm), not falling directly into the class of semilinear sectorial equations because 
the nonlinear term in it is of the same order, or value, as the main part operator. Our recent paper 
EH was devoted to such type problem, the quasi-geostrophic equation in R 2 (e.g. [10 H m [Ml 133). 
Another paper |10| was devoted to fractional Hamilton-Jacobi equation in sub-critical case. But certainly 
the most celebrated example of such critical problem is the Navier-Stokes equation in dimension two. 
We will obtain and study its solutions constructed as limits of solutions to sub-critical approximations 
when a —> 0 + , where the PA operator is replaced with its fractional power — (—PA) 1+Q ,a > 0 (P is the 
projector on the space of divergence-free functions; see e.g. EH). The 3-D N-S equation is super-critical 
in a sense specified in subsection 11.41 A higher order diffusion term like — (— PA) s ,s > |, is needed in 
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3-D to guarantee global in time existence and regularity of solutions to the corresponding regularization 
of the N-S equation. 

The classical 3 — D Navier-Stokes equation considered here has the form: 

ut = vAu — 'S/p — (u ■ S7)u + /, divu = 0, x £ fl, t > 0, 

(1.1) u = 0, t > 0, x £ <9D, 

«(0,x) = uq (x) , 

where v > 0 is the viscosity coefficient, u = (u\(t,x),U 2 (t,x),U 3 (t,x)) denotes velocity, p = p(t,x) 
pressure, and / = f 2 {x), fs(x)) external force, and Q is a bounded domain with C 2 boundary. It 

is impossible to recall even the most important results devoted to that problem, since the corresponding 
literature is too large; see anyway [14, 15j [T7J |T5J (251 [31, 32, [34} |2] together with the references cited 
there. 

In space dimensions 2 and 3 the N-S equation possess local in time regular solutions, as stated in 
Theorems ll.5l and ll.6l We analyze further criticality of the Navier-Stokes equation (compare [HUGS), i n a 
sense that available for it L 2 (Cl) a priori estimate (11.111) is not sufficient to control its nonlinearity through 
the viscosity term vPAu. Consequently, possible is a balance between the income from nonlinearity and 
the stabilizing action of the viscosity so that the local solutions can not be extended globally in time. For 
small initial data the decisive role is played by viscosity, while for larger initial data the nonlinear term 
is strong enough to destroy regularization of the solutions through the main part operator. We will try 
to see such effect through the estimates obtained in the paper. 

In 2-D, our idea is to improving a bit the viscosity term to make the whole problem subcritical , such 
that the improved viscosity together with the known L 2 (Sl) a priori estimate will control the nonlinear 
term. The way for obtaining such effect is to replace the classical viscosity term Au = —i/PAu through 
a bit higher fractional diffusion — A 1+a u, with small a > 0. Next, we will study the process of letting a 
to 0 + ; which properties/estimates of the solutions of regularized problems (with a > 0) are lost in such 
a limit. We also look at 3-D N-S equation as a super-critical problem in a sense that a stronger diffusion 
term, with —PA operator in the power s > -, is needed to guarantee the control on the nonlinear term 
with the use of the standard L 2 (D) a priori estimate. 

In 2-D, instead of m we consider a family of sub-critical problems , with a £ (0, ^]: 

ut = —A 1+a u — P(u ■ S7)u + Pf, t > 0, 

(1.2) 

u( 0, x) = u o(x). 

The approximation proposed in 3-D is given in (12.61) with s > |. 

Until very recently, fractional power operators were not seriously studied in the literature, therefore 
we hope the regularizations proposed will help to understand the difficulties faced in the original N-S 
equation. An analogous phenomena was studied recently El for the quasi-geostrophic equation in R 2 
(e.g. mEunnsniEZ])- The technique used in the present paper is similar to that of E- 

There were several tries of replacing the classical N-S equation, or the viscosity term in it, with another 
equation having better properties of solutions, starting with J. Leray a-regularization reported in paper 
[25j . see also |14l . Modification of one factor in nonlinearity, using mollifier, was sufficient to improve 
properties of solutions. Another modification of the N-S equation was proposed by J.-L. Lions in m 
Chapter 1, Remarque 6.11], where the —A operator was replaced with —A + k{— A) ; , l > |. Further 
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modifications can be easily found in the literature. In fact, we follow here the idea of J.-L. Lions to 
replace the diffusion term with a stronger one fractal diffusion term. 

1.1. Introductory facts. Notation. We are using standard notation for Sobolev spaces. Compare [33] or 
GO Chapter 1] for properties of fractional order Sobolev spaces; see also [21] for Sobolev type embeddings. 
For r£l, let r~ denotes a number strictly less than r but close to it. Similarly, r + > r and r + close 
to r. When needed for clarity of the presentation, we mark the dependence of the solution u of (E2D on 
a £ (0, i], calling it u a , or u e for the approximation in 3-D. 

Recall first [3] Chapter 3], that studying in a Banach space X an abstract Cauchy’s problem with 
sectorial positive operator A and solutions varying in the phase space = D(A^)\ 

u t + Au = F(u), 

(L3) 

u(0) = MO, 

knowing an a priori estimate of all its potential X@ solutions; ||u(t)||y < const in another Banach space 
D(A) C Y, we say that the nonlinear term F is sub-critical relative to such a priori bound , if for each 
such solution u(t) an estimate is valid 

(1.4) 

^06 (0,1 Y^nondecreasingg\ [0,oo)—>•[(), oo) \\F(umx< 9 (\\u(t)\\y)(l + \\u(t)\\%), 

for all t £ (0, t U q ), where t uo is the ’life time’ of that solution. In that case (e.g. [3] Chapter 3]) the XP 
norm of the local solution will be bounded on [0,t uo ), which allows to extend such solution globally in 
time. If the above estimate is possible for 0 = 1, but not for 9 < 1, the nonlinearity F is called critical 
relative to that a priori estimate. 

Note that critical nonlinearities are ’of the same order’ in the equation as the main part operator A 
(compare [34]). The main part operator A will not control the nonlinearity in that case, unless we find a 
better a priori estimate. 


1.2. Properties of the operator —A. Familiar in the theory of the N-S equation are the following 
spaces: 

C r (n) = [17(D)] 3 , 

(1.5) W 2 ’ r (D) = [W 2 ’ r (D)] 3 , 

x r = d £ r (n) {<() (E [C^°(D)] 3 ; div(f) = 0}, 


1 < r < oo. We define also the Stokes operator 


A r = — vP r 


A 0 0 
0 A 0 
0 0 A 


where P r denotes the projection from £ r (D) to X r given by the decomposition of /7(D) onto the space 
of divergence-free vector fields and scalar-function gradient (e.g. EE]). It is further known m Lemma 
1.1], that 


Proposition 1.1. The operator —A r considered with the domain 

D{A r ) = X r n {cj> G W 2 ’ r (D); </> = 0 on <9D}, 
generates on X r an analytic semigroup {e~ tAr } for arbitrary 1 < r < oo. 
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A complete description of the domains of fractional powers of the Stokes operator A r = —vP r A, 
D{A r ) = D(-A)nl r , can be found in [121 P- 269], or in [l 8 |. Note further, that the domains of negative 
powers of the operator A r are introduced through the relation (e.g. [191 p.269]): 

D(A>) = D(A^Y, i + i = l. 

Thus, D(A^),/3 < 0, is the completion of X r under the norm || ■ || 0 , r '- 

It is also easy to see, that the resolvent of the operator A 2 fulfills an estimate: 3 ff(er(A 2 )) > v Ai, where 
Ai is the first positive eigenvalue of —A in L 2 (B) considered with the Dirichlet boundary condition. The 
same estimate remains valid for the operator considered in L r (Q) with any r G (1, 00 ). It follows further 
from m Lem. 3.1], that the resolvent of A r is compact, also the embeddings D(A@) C IA(A“) are 
compact when 0 < a < f3 m Th. 1.4.8]). In fact the operator A = A 2 (we skip the subscript further 
for simplicity) is self-adjoint in the Hilbert space £ 2 (B); see e.g. [T7l [T9l■ For such type operators the 
powers of the order (1 + a) have similar properties; in particular they are also sectorial operators. 

Consequently, the operators A r , 1 < r < 00 , are sectorial positive.Fractional powers of the order 1 A a 
for such operators are introduced through the Balakrishnan formula (j23|): 

(1.6) A 1+a (j)= 2sm ( 7ra ) f A Q [A(A + A)- 1 ] 2 </.dA, </> G D(A 2 ), 

a7r J 0 

We recall that in case of the N-S equation the, specific for that problem, a priori estimate is obtained 
multiplying m in [T 2 (fl )] 3 through u , to obtain: 

(1.7) ^\H\cHn) = -^||Vu||^ (Q) - J^Vp-udx + J f ■ udx , 

since the nonlinear component vanish in that calculation due to condition divu = 0 : 



The term Vp • udx , for regular solutions, is transformed as follows: 


(1.9) 


Vp ■ udx 



— / pdivudx = 0 . 

J o 


Consequently, an £ 2 (B) estimate is obtained: 

( 1 . 10 ) 

\ d f IA V 

9 ~jj\\ u \\c 2 (n) = ~^ll^ M llc 2 (n) + f ' u ^ x ^ — T 7 .l^ 7 ' u llic 2 (n) + c ^ll/llc 2 (n) ^ — ^HMl£ 2 (n) + c ^ll/llc 2 (n)> 

thanks to the Poincare inequality. We finally obtain a global in time estimate of the solution: 

2c v cp\\f\\ 2 r2tri ' 

(1-11) IK*)||£ 2 (n) < max{||uo||^ 2(n) ;---— }, 

where cp denotes the constant in the Poincare inequality. Having already the last estimate one can return 
to ( 11 . 101 ) to see that 


(1-12) ll' U lli 2 (0,T;[ffi(n)] 3 ) “ 2o ( Cl/ ^ 1 |I^H'C 2 ( fi ) + H u °llic 2 (n))7 

for arbitrary T > 0. These are the strongest natural a priori estimates that can be obtained for, sufficiently 
regular, solutions of the N-S equation. 
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Remark 1.2. Note that similar estimates are also valid for solutions of the fractal approximations 
& with only one difference that the [Rg(r2 )] 3 norm is replaced with the [Rg + “(r2 )] 3 norm. Conse¬ 
quently, the estimates (II. Ill) and ( 11 . 121 ) are valid for all the solutions u a uniformly in a > 0 . 

1.3. Local in time solvability of the 3-D and 2-D N-S problems. We will rewrite N-S equation in 
a form of an abstract parabolic equation with sectorial positive operator and solve it using Dan Henry’s 
approach (M). I want to make a comment here, that using Henry’s approach we have fairly large 
choice of phase space (= space where solution varies). In fact, calling here an interesting paper [29], as 
far as we consider the corresponding to (11.31) linear Cauchy problem with sectorial positive operator: 

Ut + Au = 0, t > 0, 
u{ 0 ) = u 0 , 

we can ’set it’ at any level of the fractional power scale X® = D((— A)P),j3 £ R, corresponding to —A 
(see [l] Section V.2] for an extension of that idea). When we move to the semilinear problem (11.31) with 
nonlinearity F subordinated to A, it is an art to choose the proper level at that scale to be the phase 
space for semilinear problem. For that, we need to consider a priori estimates available for the specific 
equation, usually of physical origin, e.g. following from energy decay or conservation of mass valid in the 
process described through the equation. The full semilinear problem will be next written abstractly as 

Ut + Au = F(u), t > 0, 

■u(O) = Uq , 

The standard way to set the problem in the above setting, in L 2 (Vl) (see e.g. [31][I7[ (T5] ), is to apply 
to the equation the projector P = P 2 : [L 2 (H)| 3 —>• H , where H is the closure in [L 2 (fl)] 3 of the set 
of divergence free functions {it £ [Cq° (fl)} 3 ; divu = 0}. The pressure term disappears then from the 
equation. The realization A of the diagonal matrix operator vP[— A] 3 X 3 acts from D(A) —> H. We 
also introduce the energy space V = {u £ [H^ (ST)] 3 ; divu = 0}, and the simplified notation for the 
nonlinearity: F{u) = —P(u ■ V)u. 

Operator A = A 2 has an associated scale of fractional order spaces X& C [H 2 ^ (fl)] 3 , /3 > 0. The 
realizations of A in X 13 act from D(A' 3 ) = X@ +1 —>• X@ and are sectorial positive operators (see e.g. 

[33 H3 [201 ESI). 

We will rewrite the classical N-S equation in an equivalent form, using the property of the divergence- 
free functions. We have: 

Uj t = uAuj - + Zb dim = 0, x £ H, t > 0, j = 1,2,3, 

^ i= 1 1 

(1.13) 

u = 0 , t > 0 , x £ dfl, 
u(0,x) = u 0 (x), 
where u = (rti,H 2 ,U 3 ). 

We will recall next an estimate, important for the further calculations, borrowed from m Lemma 
2.1]. A similar observation was given also in [20] p.18] in dimension one. 

Corollary 1.3. For each j , 1 < j < N, the operator A~^P-^- extends uniquely to a bounded linear 
operator from [L r (Q)] N to X r: 1 < r < 00 . Consequently, the following estimate holds: 

II A~*P(u ■ V)u|| [L r (n) ]jv < M(r) |||u||u|||[Lr ( n)]jv. 


(1.14) 
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Observation 1.4. We have also the following estimate used further in the text. From m Lemma 
2.1] we get, for all N £ N: 


(1.15) 


||A jP(u.V)c|| [L 2 (n) p < c|||u||r>|||[ L 2 (Q)]N < c||u||[ L 4(Q)]JV ||w||[L4(0)]iV, 
\\P(u • V)v||[ L 2 (n)]jv < c||u||[ L 4( n )]w||Vv||[ L 4(n)]Jv. 


Now, for any <5 £ (0, ^), using the theory of interpolation; 


(1.16) 


\\A~ s P(u • V)i;|| [L 2 (n) p < c\\A~ip(u ■ V)w||^ 2(n)] jv \\P(u ■ V) V ||f- 2 2 ( 5 n)]N 


< c||u|| 


1+25 

[L 4 (n)r 


IIVwll 


1-25 

[i 4 (n)]N- 


In a similar way, starting from the estimates: 


\\A 2P(u-V)w|| [L 2 (Q)] iv < c|||u||v|||[ L 2 (n)]jv < c||u|| [L 6(n)]iv||i;||[ i 3(n)]jv, 

(1.17) 

\\P(u ■ V)u||[ i2(n)] N < c||u|| [i 6 (n) ]N||Vi;||[ i 3( n )]jv, 
for any 6 £ (0, we get: 


(1.18) 


\\A~ s P(u • VH| [L 2 (n)] v < c\\A~ip(u • V)u|| 2 f 2(n)r ||P(u ■ V)+||f- 2 2 ( 5 n)]N 

— C ll U ll[I/ 6 (n)] Ar H M ll[L3(f2)]iV H^ U ll[i3(f2)]AT- 


The above estimates (11.161) . (11.181) . are valid for all the space dimensions N = 2, 3,4,.... They can be 
extended further, using Sobolev type estimates, in a way depending on N. 

For local in time solvability, we will set the problem (11.131) in the base space X~i for the space 
dimension N = 2, and in the base space X~i for the space dimension N = 3. The corresponding phase 
spaces will be; Xi + C [i7 1+ (fl)] 2 ) in case N = 2, and X^ + C [P 5 + (12)] 3 ) in case N = 3 (e.g. [19], 
Proposition 1.4]). Note that, in both cases, the phase spaces are contained in the space [^“(fl)]^. Note 
also, there is another possible choice of the phase spaces (e.g. m), if we decide to work in the spaces 
[L r (Q)] N ,N = 2,3, with r > N. 

We will formulate now the corresponding local existence results for N = 2,3. 

Case N = 3. The main tool is the estimate taken from P! Lemma 2.2] (with 5 = 1,9 = p = |): 

(1.19) ||A - 8P(u • VH|[l2(q)] 3 < M||A 4 u||[ i 2( n )]3||A 4 u||[L2(n)]3- 


Since the form above is bi-linear, we have also the following consequences of the last estimate: 


( 1 . 20 ) 


||A »P((ll - v) ■ V)u||[i2(Q)]3 < M||A 4 (u — u) || [Z,2(fi)]3 ||A 4 u||[ i2 (Q)]3, 


3 + 


, 2 + . 


||A *P(m • V)(u - u)||[i2(n)]3 < M\\A* u||[ L 2( n )]3||2l 4 (-tt — -u)||[z,2(f2)]3. 


Consequently, the nonlinear term F(u) + Pf = —P{u ■ V)u + Pf acts from D(A* ) C [H% (12)] 3 into 
D(A~s) as a map, Lipschitz continuous on bounded subsets of D(A 4 ). According to [20][3], this suffices 
to obtain a local in time solution of the 3-D equation (11.11) (also of 3-D equation (12.61) 1. more precisely: 

_ i 34- 

Theorem 1.5. When Pf £ D(A »), u o £ D(Af ), then there exists a unique local in time mild 
solution u(t) to (11.11) in the phase space D(A% + ) C [P5 + (J2)] 3 . Moreover, 


( 1 . 21 ) 


u£C([0,r);P(A3 )) n C7((0, r); D(A *)), u t £ C((0, r); D(A*~)). 


Here r > 0 is the ’life time ’ of that local in time solution. Moreover, the Cauchy formula is satisfied: 


u(t) = e At uo+ ( e A ^ ^F(u(s))ds, t£[0,r), 
Jo 
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where e~ At denotes the linear semigroup corresponding to the operator A. We need also to mention that 
the considered here mild solutions have additional regularity properties, as described in particular in [3] 
p.218]/ here u G C 1 ((0, r); D(A$ )). This property is used in the calculations below. 

Case N = 2. We will use a version of the estimate in [T9l Lemma 2.2] (with 6 = j, 6 = p = |): 

(1.22) \\A~^P(u ■ V)u||[i 2 (o )]2 < M\\A* 'u||[l2(q>] 3||A 5 u||[L 2 (n)] 2 . 

Since the form above is bi-linear, we have also the following consequences of the last estimate: 

\\A~*P((u - v) ■ V)i;||[ L 2 (q )]2 < M\\Az (u - 'r)||[L 2 (fi)] 2 \\A* 

\\A~*P(u ■ V)(u - u)||[L 2( f2 )]2 < M\\A? (u - w)||[ i 2(n)]2||242 ( u - LOII[L 2 (n)] 2 - 


(1.23) 


Consequently, the nonlinear term F(u) + Pf = — P{u ■ V)u + Pf acts from D{A? ) C [iT 1+ (f2)] 2 into 
D(A~i) as a map, Lipschitz continuous on bounded subsets of D(A 1+ ). According to [SEEGO, this suffices 
to obtain a local in time solution of the 2-D equation (11.11) (also of 2-D equation (11.21) 1. more precisely: 

Theorem 1.6. When Pf € D(A~i), Uq € D(A^ + ) C [H 1+ (fl)] 2 , then there exists a unique local in 
time mild solution u(t) to (11.11) in the phase space D(Ai + ) C (ff 1+ (fl)] 2 . Moreover, 

(1.24) U eC([0,T);D(#))nC((0,T);D(Af)), u t G C((0,r); D{A*~)). 

Here t > 0 is the ’life time ’ of that local in time solution. Moreover, the Cauchy formula is satisfied: 


i(t) = e At uo+ ( e A ^‘ s )F(u(s))ds, tG[0,r), 
Jo 


where e At denotes the linear semigroup corresponding to the operator A. 

1.4. Super-criticality of the N-S equation in 3-D, and criticality in 2-D. We will consider now 
criticality of the N-S equation m in the sense SUE], that means ’the nonlinear term is of the same 
order as the main part operator’. 

Case N = 3. The 3-D N-S equation is super-critical in the sense stated below. A general problem 
will be called super-critical with respect to the a priori estimate in Y, if for all its possible X 13 
solutions an estimate holds 


(1.25) 


A#q>iV 0 >#q i^nondecr easing g:[0 ,oc) —>[0 ,co) \\F{u(t)\\ x <g(\\u(t)\\ Y )(l + \\A^um e x ), 


for all t G (0,t U o ), and such type estimate is not true for exponents 6 < 6q. 

In case of the Navier-Stokes equation m the mentioned above a priori estimate will be the Y = 
[L 2 (fl)] Ar estimate (11.111) . 

We will see that, as a consequence of (11.181) . the problem (11.11) is super-critical. Indeed, the estimate 
(11.181) written for the local solution u = u(t) obtained in Theorem II. 51 extends for N = 3, with the use of 
the Nirenberg-Gagliardo type estimates, to: 

(1.26) ||A 8 P(u ■ V)u||[ L 2(Q)]3 < c||u||[i6(Q)]3||u||^3^ n ^ 3 ||Vu||^3( n ^3 < C ll' U ll^ f |^ n jj 3 ll U ll[L 2 (n)]3. 

This shows the nonlinearity of the 3-D N-S equation is super-critical as a map from X~s to X~^. Exponent 
9 obtained above equals I. 
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Using the estimate of Corollarv ll.3l we will find now the value 6 0 in case of the 3-D N-S equation. We 
obtain, 

\\P(u ■ V)u|| x _i = || A~ 2 P(u • V)u||[ I/ 2 (n)] 3 < M|||u| 2 || [i 2 (n)] 3 < c||u|| 2 L 4( n )]3 


(1.27) 




and the estimates are sharp. Therefore, 9o = | in that case. 

Finally, for the further use, we will check how large the exponent s > 1 should be, for the ’strengthen 
diffusion’ of the form (—A) s replacing the usual (—A) operator in (11.11) to make the L 2 (fl) estimate 
critical. Using again Corollary 11.31 we get that 

||A 2 P(u ■ V)ii||[£ 2 (Q )]3 < M\\\u\ |[L 2 (fi )] 3 A cIMIjif 23 "-^^] 3 ll u ll[L 2 (n)] 3- 


(1.28) 


We find such critical value of s from the condition; 2 { 2 s-\) = -*-■ Consequently, s = |. 

Case N = 2. Using Corollary II.31 it is easy to see criticality of the 2-D Navier-Stokes equation. More 
precisely, to verify that its nonlinearity is critical with respect to the standard L 2 (fl) a priori estimates, 
as a map from Xi C [i7 1 (U )] 2 to X~i. Indeed, 


(1.29) 


||P(u-V)u|| i < ||A 2 P(u • V)u|| [i 2 (n)]2 < M|||u| 2 || [L 2 (n)] 2 


< C IMI[Z>(fi)] 2 < l|w||[_f/l(n)] 2 l|M||[L 2( n )]2 = c(||u||[ L 2 ( n )] 2 )||u||[ J? l (n)] 2 , 
and no better estimate (with exponent smaller than 1) is possible. 

Observation 1.7. Using the estimate (11.161) with S = j, we will show that the nonlinearity in 2-D 
N-S equation is critical as a map from Xi C [Hi (U )] 2 to X~i. Indeed, from (11.161) with S = j, N = 2, 

(1.30) ||A 4 P( m . V)u|| [L 2 (n)]2 < c||u||[ 2 z , 4 (Q)] 2 ||Vu || [ 2 L4(n)]2 < c||u|| [L 2 (n)] 2 ||u|| [ff 3 (n)]2 , 

where the Nirenberg-Gagliardo type estimates were used: 


(1.31) 

No better estimates are possible. 


II0IIz. 4 (O) < c ll^ll^§ (a) H < ? !, ll l 2 ( n)> 
I0llw 1 .4(n) < c|| 0 || s 


h ?( n) 


Note further, that in 2-D the nonlinearity is also critical with respect to the L 2 (U) a priori estimates 
as a map between X 1 = D(A) C [H 2 (H)} 2 and X 2 . Indeed, the following estimate holds, 

(!- 32 ) ||P(u • V)u||x 2 < c||u||[ L 4 (n) ]2||Vu|| [i 4 (n)] 2 < c||u||[ ff 2 (n)] 2 ||u||[ i 2 (n)] 2 , 

and no smaller exponent on the H 2 {fX) norm is possible. 


2. Global in time solutions in 3-D. Small data. 

As well known, global in time extendibility of the local mild solution constructed in the Theorem 11.51 
is possible provided we have sufficiently well a priori estimates that prevents the D(Ai + ) C (ff5 + (f2 )] 3 
norm of the solution to blow up in a finite time. We will propose next such type estimate, in [Af 1 (U2)] 3 , 
for solutions of the 3-D N-S equation when the data; uq and / are smooth and sufficiently small. Such 
estimate will be used later to construct global in time solutions of the N-S with small data. Another 
approach to that problem in arbitrary dimension N , using estimates on integral equation, was presented 
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Theorem 2.1. If uq £ D{A) C [TT 2 (ST)] 3 and f £ [L 2 (n)] 3 fulfill the ’smallness restriction’ (12.41) . 
then the [P 1 )!!)] 3 norms of the solutions u are bounded uniformly in time t > 0. 


Proof. It is known (e.g. [Uli), that the local in time solutions u are regularized for t > 0 (through the 
equation) into the space D(A) C [P 2 (D)] 3 . When the initial data uq £ D(A ), they simply vary in D(A ) 
for t > 0 small, until possible blow-up time t(a, uq). We want to show that, to small uq and / correspond 
small solutions, in the [P^fl )] 3 norm (uniformly in time). 

To get estimate of the solution u in [I? 1 (fi )] 3 multiply (11.21) by Au to obtain: 


(2.1) < u t ,Au >[z,2(n)]3= - < Au,Au >[z,2(n)]3 - < P(u-V)u,Au >[l 2 (o)] 3 + < Pf,Au >[L2(n)]3, 

since the pressure term vanishes. Thanks to (11.151) . this gives, 


1 d 


9 ^IMI[Hl(n )]3 < - C HMI[H 2 (n) ]3 +c||u||[WL4 ( n)]3||u|| [L 4 (n)] 3||u|| [ff 2 (n)] 3 + ||P/||[L2 (n)] 3||u|| [ff 2 (n)] 3 


( 2 . 2 ) 


< —cv ||u||fi/ 2 (n)]3 + c i|l u ll[ 4 fl-i(n)]3ll' u ll[ J H-2( n )]3 + \\P f\\[L 2 (n)] 3 IMI[.H 2 (n)]3 


< - 


C V . 


2 iMIf-spo )] 3 + ^(ll M ll[l/i(n )] 3 + \\Pf\\[L 2 (n)] 3 )> 

with the standard use of Young’s inequality, and the embedding H 2 (ffL) C H 1 ( fl) (constant c'). Denoting; 
y(t) := || w (0llfffi(n)] 3 ) we arrive at the differential inequality (e.g. [30l 135] ): 

^ ^ + C v y*{t) + C' l/ ||P/|| 2 i 2 (Q)] 3 , 

y{ o) = ll M (o)llfi/i(n)] 3 - 

Analyzing its right hand side, real function g(z) = — ^ffz + C v z 3 + C V \\P f\\ 2 L2 ^^ 3 , we see that g( 0) = 

/ 1 

^‘•'\\Pf\\[L 2 (n)] 3 > < 0, and g has a minimum for the argument z min = (j§£) 4 , with g(z min ) < 0 

when the ’free term’ C' !/ ||P/||? L 2 (qv 3 is small. More precisely, to keep the value of y(t) bounded for all 
positive times, we need to assume the smallness hypothesis: Let the data: ||tt(0)|| and ||P/||f L 2 (n)] 3 

be so small, that: 


(2.4) 


Q^min ) < 0, equivalently \\P f\\ [L 2 {n)]3 < 


.. ||2 f 3 c'v 

ll U 0||[JJl(n)] 3 < Zmin = I 


2 2cV 

^ and 


7CJ 


Note that 4 is proportional to v 2 . Consequently we obtain the bound 
(2.5) IH^IIfi/qn )] 3 < z min valid for all t > 0. 

With the last assumption, the smooth local solutions u(t), introduced in Theorem 13.II are bounded in 
[P 1 (D )] 3 uniformly in t > 0. Note, that a bit more accurate bounds in (12.41) are possible if one compare 
the data with the two positive zeros of the function g. See also the corresponding restrictions formulated 
in [311 Theorem 3.7]. □ 


2.1. Regularization of the 3-D N-S equation. In the 3-D case the N-S equation ED is super¬ 
critical. Global in time extendibility of the local solutions constructed in Theorem 11.51 is in general not 
possible (unless for small data). The viscosity term in the equation (11.11) is not strong enough to control 
the nonlinearity and the regularization used in 2-D case is also not sufficient to work in the 3-D case. In 
the present subsection we propose another approximation/regularization (with stronger diffusion term) 
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of the 3-D N-S equation, for which the L 2 (D) estimates (same as for the original equation) are sufficient 
to make such problems subcritical. 

We propose now approximation of the original 3-D Navier-Stokes equation having global in time, 
unique and regular solutions. Consider namely the approximation/regularization of (11.11) of the form: 


( 2 . 6 ) 


ut = — (A H——j-A s )w — P(u ■ V)u + Pf , t > 0, 
w(0, x) = u 0 (x), 


with a parameter s > 1 (to be chosen), and e > 0. For fixed (for a moment) parameter e, denote the 
solution to the above problem as u s . We mean here the solution on the base space D(A~i), obtained in 
a similar way as for the original N-S equation (ED in subsection [Q1 

It is clear how to determine the proper, sufficiently large, value of the exponent s > 1 to guarantee, 
together with the L 2 (D) estimates valid also for solutions of (12.61) . that the nonlinear term is subordinated 
to the main part operator (A + -^f=tA s ). We need to compare the bound obtained for the nonlinearity 
with the income from the improved viscosity term. Estimate of the nonlinearity, obtained from (11.161) . 
reads 


( 2 - 7 ) II- 4 * F {u s {t)) ||[i2 (n )]3 < c||M s (t)|| [ 2 i4(n)] 3||M s (t)||[V 1 , 4( n )] 3- 

Let s > |. With the use of the Nirenberg-Gagliardo type estimate we obtain the subordination condition 


( 2 . 8 ) 

|| J 4 _I F(u s (t))||[ i 2(n)]3 < c(||w s (t)||[ i 2 (n) ]3)|| u s (t)||^ 2s _i (n)]3 


< C '(||u S (t)|| [i 2 (n)] 3)||^ 


[L 2 (n)]3, 


where 0 = < 1- Consequently, the standard L 2 (D) a priori estimate is sufficient to assure the global 

in time extendibility of such local solutions u s when s > |. We will study such approximation next. 


Remark 2.2. The exponent ^ proposed for regularization of m in m Chapter 1, Remarque 6.11] 
is the same (also, as in (11.281) ). 


To make further considerations simpler, let us consider the introduced above regularization of the 
original 3-D N-S equation with exponent s = 2 in (12.61) (which is a particular choice). We will describe 
next shortly the process of passing to the limit, as e —> 0 + , in such approximations. The idea is similar as 
in nznm using the parabolic regularization technique. An abstract counterpart of such regularization 
(s = 2) of the 3-D N-S equation has the form: 

(2.9) ut = — Au - A 2 u — P(u ■ V)w + Pf, it(0) = uq. 

Solutions of that problem will be denoted further as u e . 

We formulate next the corresponding local existence result for such problem. 

Theorem 2.3. Consider (12.61) with s = 2 as equation in D(A~i). When Pf £ D(A~i), uq £ 
D(A%) c [if 3 (fl)] 3 , then there exists a unique local in time mild solution u e (t) to (12.6|1 . s = 2, in the 
phase space D(A%). Moreover, 

(2.10) u £ eC([0,r);T(Af))nC((0,T);T(Ai)), u\ £ <7((0, r); D(A*~)). 

Here t > 0 is the ’life time’ of that local in time solution. Also, the corresponding Cauchy formula is 
satisfied. 
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Moreover, together with the standard [L 2 ( f2)] 3 a priori estimate valid for solutions of (12.61) uniformly 
in e > 0, the obtained above local solutions will be extended globally in time in the class (12.101) . 


The proof, similar as presented previously, is omitted. We only resolve the L 2 estimate for (12.91) . with 
an extra term compare to (11.101) . Multiplying (12.91) by it e we obtain 

< ul,u e >[i,2(Q)]3= - < Au e ,u e >[z,2( n )]3 -- < A 2 u e ,u e >[£2( 0 )]3 

( 2 . 11 ) 


- < P{lT ■ V)lt e ,U £ >[L2(Q)]3 + < Pf,U e >[L2( 0 )]3 . 


Noting that the nonlinear term vanishes, we get an estimate 

(2-12) - ^||M e (i)||jL2(Q)]3 < -\\A 2 U <i \\ 2 L 2^3 — -||^M £ ||^ L 2(n)]3 + ||^ > /||[L2(n)]3||M e ||[L2(n)]3- 

The reasoning goes as for (11.101) . if we neglect the second right hand side term, giving precisely (11.111) for 
u e . In addition to (11.121) . for arbitrary fixed T > 0 we obtain an extra estimate 

(2-13) e\\u' |li 2 (0,T;[H 2 (O)] 3 ) ^ 2 + 

The latter estimate will be used to pass to the limit, when e -A 0 + , in the weak formulation of (12.91) . The 
idea is taken from [27[ Chapter 3]. 

The problem (12.91) has a weak formulation; for arbitrary test function v £ D(A) and fixed arbitrary 
T > 0, we set 


(2.14) 


d 

-T<U,V > [Z ,2 (n)] 3 = 


< A^u,A^v >[L2(n)]3 - - < Au,Av >[£2(n)]3 

- < P(u ■ V)u,v >[£2(q)]3 + < Pf,v >[L 2 (n)]3, 


and the weak solutions are expected in L 2 (0,T; D(A)). The solutions u e constructed in Theorem [2T3] fulfill 
the above weak form of dZSD - The time derivative here will be understand in the sense of the ’scalar 
distributions’ (e.g. Ea E]). Passing to the limit in the nonlinear term is possible thanks to the uniform 
estimates of u e in L 2 (0, T ; [Hq (Id)] 3 ), as described below in 2-D case. Note also that, when letting e —> 0 + , 
the second right hand side term will tend to zero in the space L 2 (0,T; [-/^(fl)] 3 ), thanks to (12.131) . An 
even more interesting observation is that, when passed to such a limit, we obtain a ’weak solutions’ of 
the original N-S equation (e.g. m Chapter 3.1]), due to J. Leray. The only difference is that the set of 
test functions v £ D(A) we are using is smaller. But it is dense in D(A^), so the two formulations are 
equivalent. We will not extend that considerations here in details. 


3. 2-D CRITICAL N-S EQUATION AS A LIMIT OF SUB-CRITICAL APPROXIMATIONS. 

We will describe now the convergence of the solutions of the fractal approximations (11.21) to the solution 
of the limiting 2-D N-S equation. Precisely as in Theorem 11.61 the solutions u a of (11.21) (the superscript 
is added for clarity) will be constructed in the class: 

(3.1) u a £ C([0,ry,D(Ai + ))nC((0,ry,D(Ai)), < £ C((0, r); D(Af)). 

Further, the L 2 (fl) a priori estimates are satisfied for u a uniformly in a £ (0, |]. Consequently, we claim: 

Theorem 3.1. The local in time solution u a (t) of the approximating problem (11.21) constructed as in 
Theorem \m will be extended globally in time in the above class. Moreover, the standard L 2 (fl) a priori 
estimate (11.111) is valid for it uniformly in a £ (0, ^]. 
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All the nice properties of the abstract semilinear sectorial equation, like elegant theory of existence 
and uniqueness, regularization of the solution for positive times, are valid for solutions u a (t) with any 
a > 0. Typical regularization goes from the phase space X^,fd < 1, to the space X 1 = D(A); see e.g. 
[20113]. However, for problems m, estimates of norms better that L 2 (fl) will depend on a > 0, possibly 
blowing up when a —> 0 + , and therefore can not be extended in general to a limit solution u. 

3.1. 2-D (ED as a limit of (11.21) when a —>• 0 + . A precise description of letting a —> 0 + in the 
equation m is given next. In this section we consider the solutions u a of m constructed in Theorem 
IQ on the phase space D(A *) C [H 2 (H)] 2 . Such solutions, for any a £ (0, |], are varying in [H 2 (H)] 2 . 
Indeed, according to m Proposition 1.4]: For any /3 > 0, the domain D^A^) is continuously embedded 
in X 2 D [H 2 P(VL)\ 2 . Solutions u a fulfill also, uniformly in a £ (0, |], estimate (11.111) in [L 2 (Q)] 2 . More 
precisely, for such solutions of m we have an estimate: 

(3-2) ^corasi>0^ae(0,i] ll M I|l° 0 ([0,oo);[L 2 (0)] 2 ) A COUSt. 

This is the main information allowing us to let a —> 0 + in the equation US- 

Passing in (11.21) to the limit a —> 0 + . We look at OD as an equation in [L 2 (fl)] 2 , and ’multiply’ by 
the test function A~ 1 ~ a (j) where <f> £ D(A%); recall that D(A%) C [I?5(H)] 2 , 

(3.3) 

<v? + P(u a ■ V)«", A- 1 -*# > [iW = - < A 1+a u a , A _1_Q </> > [LW + < Pf, > [£W . 

We will discuss now the convergence of the terms in (13.31) one by one. Note that when a —> 0 + 
then, by Lemma 14.11 A _1_a <)> —> A _1 </>. Thanks to uniform in a £ (0, A boundedness of u a in 
L°°([0, 00 ); [L 2 (H)] 2 ), we obtain: 

< A 1+a U a , A _1_a <^ >[L 2 (n) ]2 = < u a , (f> >[L 2 (fi)] 2 ^< M, (j> >[L 2 (n)] 2 , 

< Pf,A~ 1 ~ a (j) >[L 2 (n)]2-K >[L2( n )]2, 

where u is the weak limit of u a in [L 2 (fl)] 2 as a —> 0 + (over a sequence {a n } convergent to 0 + ; various 
sequences may lead to various weak limits). 

We return to (El to see that letting a —>• 0 + over a sequence {a n }, where u denotes weak limit in 
[L 2 (H)] 2 of such sequence, we have 

(3.5) lim < A _Q (u“ + P(u a ■ V)n Q ), A^ 1 ^ >[i, 2 (n)p = — <u,(j> >[L 2 (n)] 2 + < Pf^A^ 1 ^ >[z, 2 (Q)p, 
since the right hand side is convergent. Consequently, the left hand side has a limit as a n —> 0 + , 

(3.6) lim < A”“(u“ + P(u“ • V)u“), >[L 2 (n)] 2 = 

Ot-n ^0 

Note that, for the ’test functions’ A^cf) varying in a separable Banach space, passing countable many 
times to a subsequence, we can chose a common subsequence proper for all test functions in dense subset 
of the space. Consequently, the equation below will be fulfilled in the whole space. We obtain 

(3-7) ^<t> = - < U, (j> >[L2(f2)]2 + < Pf, A~ l (j) >[L2(Q)]2, 

which is a weak form of the limiting equation. 

Separation of terms. The two terms of [A _a (u“ + P{u a ■ V)u“)] will be separated when letting 
a —> 0 + . More precisely we have: 
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Remark 3.2. Since the approximating solutions u a satisfy (in particular) 
(3.8) u a G L°°(0, T\ [L\n)} 2 ) : < G L\ 0, T; [L 2 (n)} 2 ), 


then by EH Lemma 1.1, Chapt.III] 

(3-9) V^ e x 2 < Ut,r] >[l 2 ( o)] 2 = — <u a ,r] >[z, 2 (si)] 2 -> — <u 7 r] >[z. 2 (r 2 )] 2 , 

the time derivative and the convergence are understood in T>'(0,T 7 [L 2 (il)] 2 ) (space of the ’scalar 
distributions ’mv- Consequently, 


(3.10) 


= — < .4 a u,A 1 cj) >[l 2 ( o )]2 +u\, 


where is a limit in V (0,T;[L 2 (LI)} 2 ) of < P(u a ■ V)w“,A 1 a (j) >[x, 2 (fi)] 2 over a chosen sequence 
ot n -t- 0+. 

Convergence of the nonlinear term F(u a ) = P(u a • V)zt“ will be discussed next. As seen from (11.121) . 
the approximating solutions u a are bounded in L 2 (0,T; [iJ^Jl)] 2 ) uniformly in a > 0. Estimate (11.151) .' 

(3.11) ll^ -3 -F 1 (u a )||[ l 2 ( q )]2 < c||M“|| 2 i 4( n )]2 < c , ||u a ||[ L 2 (n) ] 2 ||u a ||[ ff i(o)]2,7V = 2, 


together with the consequence of equation El 




-A^u a + A~^ +a) F(u a ) + A~^ +a) Pf G L 2 (0, T; ]L 2 (Li)] 2 ), 


show that uf are bounded in L 2 (0,T] D(A x )) uniformly in a & (0, ^]. By Lions compactness lemma [271 
Theorem 5.1], the family {'u a } Q ,e(o ) ±] as bounded in the space 

w = £ L 2 (0,T; [Hg(n)] 2 ),(j) t G L 2 (0,T-,D(A~ 1 ))}, 


is precompact in L 2 (0 ,T-,[Hq (fl)] 2 ). In particular any sequence {u a "},a„ — > 0 + , has a subsequence 
convergent almost everywhere in (0 ,T) x H. Moreover, the above compactness allows to pass to a limit 
(using estimates of Observation o in the nonlinear term in (ED. Indeed, if u an —>■ u,a n —> 0 + in the 
above sense then, using a consequence of ua Lemma 2.2]; for small e > 0 

\\A-d+^P[(u • V)u - (u““ • V)u“»]|| [L 2 (n)] 2 < C\\\u - |(H + |« a »|)|| [£ . (n) ]2, 


where ^ ^ + e (so that s < 2), we have 


(3.12) 

[ < P{u a » • V)u an ,A~ 1 ~° tn (j) > [L 2 (n)] 2 dt = 

J o 


[ < A-^+^P(u an • V)u a ",A _ i-“"+V > 
Jo 

[ < A~ { i +e) P(u- V)u,A~i +t (l)> 

J 0 

The construction presented above allows us to formulate the following theorem: 


[i 2 (n)]2 dt 
[i 2 (fi)] 2 dt. 


Theorem 3.3. Let {M“} ag (o i] be the set of regular D(A%) solutions to sub-critical equations (11.21) . 
Such solutions are bounded in the space [L 2 (J2)] , uniformly in a. As a consequence of that and the 
regularity properties of such solutions (varying in D(A*) C [H^(M)] 2 ), for arbitrary sequence {a n } C 
(0, |] convergent to 0 + we can find a subsequence {«„J that the corresponding sequence of solutions 
{u an k } converges weakly in ]L 2 (fl)] 2 to a function u fulfilling the equation: 

(3-13) V<j><EX 2 Wf = — < u,<f> >[l 2 (q )]2 + < P/,A“V >[L 2 (fi)] 2 ■ 
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(3.15) 


Due to denseness of the set D(A*) in X 2 , the right hand side of (13.131) defines a unique element in X 2 . 
The left hand side is defined in (13.61) and discussed in Remark \3.2[ 

Remark 3.4. As well known (e.g. [271 Theorem 6.2]), the L°°(0, T; [L 2 (fl)] 2 ) D L 2 (0, T; D(A?)) 
solution of the 2-D N-S equation having ut £ L 2 (0,T; D(A~%)) is unique. Indeed, if u\,u 2 are two such 
(weak) solutions, applying projector P, taking the difference of the equations and multiplying the result 
by w = u\ — u 2 , we get: 

( 3 - 14 ) ^\\w(t)\\ 2 [L 2 m 2 < -\\A?w\ |f L 2 (n)] 2 + < P(w • V)ui, w > [ / / 2 ( q )] 2 + < P(u 2 • X)w,w >[l 2 ( 0 )] 2 • 

The last term vanishes for divergence-free functions. The earlier term, using an equivalent form of the 
nonlinearity in (11.131) and the Nirenberg-Gagliardo inequality, is estimated as follows: 

I < P(w ■ V)ui,W >[£ 2 (f 2)]2 I < c||w|| 2 i 4(Q)]2||A 2 Ul||[ L 2( n) ]2 

< c|k||[L2 ( Q)]2||24 2 M;||[ i 2 (n)] 2||245'Ul||[ i 2 (n)] 2. 

Inserting the last estimate into (13.141) . using Cauchy’s inequality, we obtain 

2 ffi ll w (0llfL 2 (n)] 2 < -\\ A ^ w \\[l^(q.)y + c ll w ll[L 2 (n)] 2 m 5 w||[ L 2 (n )]2||24 2 Mi||[ L 2( n) ]2 

< C|Mljx 2 (n)] 2 ll^ 2M illfL2(o)]2- 

Since ||w(0)|| 2 i2 (Q)]2 = 0 then ||w’(0llfL 2 (n)] 2 = 0 for all t £ [O,? 1 ], due to the Gronwall lemma. 

Remark 3.5. In 2-D, convergence of the approximated solutions u a of (El to the unique solution u 
of the N-S equation holds in a better sense when uq £ D(Ai). Multiplying (11.21) by Au a , we get 
(3.16) 

ua \\[ L2 (Q)] 2 = H|A +2M “llfL 2 (fi)] 2 + < P(u a ■ X)u a ,Au a >[z, 2 (q )]2 + < Pf,Au a >[z, 2 (n )]2 . 
Further, due to (11.161) . Nirenberg-Gagliardo and Young inequalities, 

| < P[u a ■ V)u a ,Au a >[£3(n)]2 | < \\A~^P(u a • V)u Q ||[L 2 (n )]2 H^ 1 " 1 " 2 u“||[£ 2 (q )]2 

2[l+g-2a 2 ] 2[2 + 3Q- + a 2 ] 

— c ll uQ ll[’zI" 4< (n)] 2 ll^ ua llfL 4 '(n)] 2 ll^ 1+2MQ; ll['£' 2 (^)] 2 — e ll^ 1+2MQ llfL 2 (n)] 2 + C e ||u“||[ ff i 1 (Q®“ 2 l|w“ll[ L 2(Q)] a 2 • 

For small e > 0, a > 0, we thus have 

1 J I 2[l+n-2a 2 ] 2[2 + 3f + a: 2 ] 

( 3 - 17 ) 2 dt ll 242uQ H? i2 ( n )] 2 - ~^W Al+2 u “llfi 2 (n)] 2 + c, e||w a ||[ i/ il (^ 5 ) “ 2 KH[^ + ll-P/II^L 2 (n)] 2 - 

It follows from the natural a priori estimates that the family {it Q } a g(o i], is bounded in the norm 

(3.18) L°°(0, T; [H^D)] 2 ) n T 2 (0, T; D(A)). 

Consequently, it is precompact in 

(3.19) 2/(0, T; [R 0 1_ (D)] 2 ) D L 2 (f), T; D{A 1 ~)), 

with arbitrary p > 1. The convergence u a —> u as a —> 0 + is thus verified in the above space. 


Remark 3.6. To understand better the utility of the method used in the study of the global extendibility 


of solutions we will discuss shortly an 




of the Burgers type system in 3-D, obtained by neglecting 


^suggested kindly by professor Eduard Feireisl 
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the viscosity in homogeneous 3-D N-S equation: 

U t = vAU - (U • V)Z7, xGCt,t> 0, 

(3.20) U = 0, t> 0, xedfl, 

U(0,x) = U o (x). 

It is easy to see that each component of the sufficiently regular (that is, varying in [L°°(f2 )] 3 ) solution of 
(13.201) fulfills Maximum Principle: 


(3.22) 


(3-21) \\Ui(t, •)IU~(Q) < ||CA)i||L«(fi)) i — 1,2,3. 

We are thus given a natural a priori estimate inY= [L°°(r2)] 3 for such system. The nonlinear term is 
as in the N-S equation, and we have the estimates: 

\\{U ■ V)C/||[ L 2(n)]3 < ||C/|| [i°°(f2)]3 ||C/|| [TJ-1(0)]3 < ||t>o|| [i°°(f2)]3 l|C/||[TTi(f2)]3 5 
II • v)c/’|| [ff -i (n) ]3 < c||tr||f ioo(n)] 3, 
valid in particular for all the local solutions varying in [R5 + (D)] 3 fl [i/^fi)] 3 (as in Theorem !!. 5\) . Since 

3 d - _____ 

H 2 (D) c the Maximum Principle works for such solutions. By interpolation, (13.221) gives 

(3.23) ||(C/-V)C/||^_ 4(f2)]3 < cllC/ll^^jallC/ll^^a < c||C/ 0 ||^ (£7)] 3||C/||^ ol(r2)] 3, 

which shows the nonlinearity is sub-critical in that case. Consequently, the local solution will be extended 
globally in time. 

The above example shows that our approach is sensitive not only on the form of nonlinearity, but also 
another specific properties. It also indicates the role of the pressure in the classical N-S equation, the 
term which seems responsible for the delicate properties of solutions of that equation. 


Remark 3.7. It is evident from the considerations above that the phenomenon of loosing regularity by 
local smooth solutions of 3-D N-S is possible only if they enter the super-critical range of the nonlinearity 
(compare (11.261) ). It would be therefore interesting to consider local solutions corresponding to initial data 
u o fulfilling 

(3-24) ||y4^8F(M 0 )||[L2(O)]3 > g(||Uo||[ L 2( n) ]3)||A8M 0 ||^ h 2 f ( 0 )]3, 

where 0 < e < as eventual candidates for such phenomenon. 

4. Appendix. Fractional powers operators and estimates. 


4.1. Some technicalities. When passing to the limit in the considerations above it was important that 
the estimates, in particular the constants in it, can be taken uniform in a. Therefore, in the technical 
lemmas below we need to care on a very precise expression of that uniformity. Even some estimates can 
be found in the literature, usually such uniformity is not clear from the presentation, thus we include 
here the proofs for completeness. 

First, we formulate a lemma used in the previous calculations: 

Lemma 4.1. Let A be a positive operator in a Banach space X nsaiMiG]). For arbitrary (f> £ X , 
we have 

V £>0 3 l ||(/ - A-^cfWx < sin(7r/j)( 2I/(1 + M) + ||. Y + e. 

7r 

Consequently, the left hand side tends to zero as 0 < /3 —» 0 + . 
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Proof. Our task is, for fixed (f> £ X and /? near 0 + , to estimate the expression: 


{A-? - I)cf = 


sin(7 t/3) 


A _/3 (A + A)~ 1 (/)dX - 


sin(7r(l — (3)) f°° A^ 1 ^ 1 


(4.1) 


A + 1 


dX (j> 


f°° A“^[(A + - -l—cfldX. 

TT JO A + 1 

/ing 

f 

Jo 


In the estimates we are using the following properties; taken from [2H1 P- 62] equality valid for rj £ (0,1) 

roo ^-1 


-dX = 


the simple formula: 
(4.2) 


A + 1 singly) 


1 


(A + A) 1 (f — . </> — . 1 [A(A + A) 1 f> - cj) + (A + A) 


-i. 


x + r a + i 

and the two asymptotic properties of non-negative operators valid on functions cj) £ X taken from 
Proposition 1.1.3]: 


(4.3) 


lim A(A + A) 1 cf> = ((>, 

X—too 

lim (A + A)~ 1 A(f> = 0. 

A—>-oo 


Returning to the proof, we split the integral in (14.11) into (0, L) and (L, oo) and estimate the first part, 


sin(7T/3) 


(4.4) 


[ L A-^(-^- - (A + Ay^dXWx < [ L X~P(1 + M)dX\\(j)\\x 

Jo A + 1 TT Jo 


sin(7r/3) L 1 ^ 


(1 + M)U\\ X , 


7T 1 - , 

where L > 0 will be chosen later. Note that letting j3 —> 0 + the result of the estimate above is bounded 
by | Sln ^ r<3 ) |2L(1 + M)||^>||x —> 0, for any fixed L > 0. 

Next using m , the integral over (L, oo) is, for <f> £ X, estimated as follows: 

TT J L A + l tt X + 3 

where due to (14.31) we see that 


-||A(A + A) 1 4>— <f>+ (A + A) 1 (f>\\xdX, 


M 


(4.6) || (A + 1)(A + A) A<(> — </ , ||x<||A(A + A) (f> — </>||x + || (A + A) </>||x is £ + z —;—r ||<^|| x as A —> oo, 


1 + A 


e > 0 arbitrary fixed. Consequently we obtain: 

(4.7) 


sin(? t/ 3) f°° X 0 , sin(7r/3) 

- / TTrl 6 + Tyr r? A' “A < - 

TT J L X + l' 1 + A 7T 


M , 


A P 1 (e + — H0Hoc)rfA 


< + iAAAXAuuu, 


7T p 7T /3 + 1 

for sufficiently large value of L > 1, as specified in (14.61) . Note that letting fj —> 0 + in the resulting 
estimate we have: 


(4.8) 


sin(7 t 8) r _ R sin(7r/3) L 13 1 . 

- 4-^-L 0 e+ ---—— M 0 \\ x < e + sm(n/3)L M 0 \\ x , 

7T p 7T P + 1 


for chosen large value of L. 

For such L we get a final estimate of the integral in (14.11) having the form: 

(4J) 

WiA-P-mx < ^^^11 r X-^X+Ar^-j^-^dMlx < S m(7rP) C L{1 + M) +L 

TT J Q A + 1 TT 


x+e, 
































NEW LOOK AT THE NAVIER-STOKES EQUATION 


17 


where e > 0 was arbitrary. The right hand side of (14.91) will be made small when we let fi near 0 + , noting 
e was an arbitrary positive number. □ 


4.2. Properties of the fractional powers operators. Recall first the Balakrishnan definition of frac¬ 
tional power of non-negative operator (e.g. [23j p. 299]). Let A be a closed linear densely defined operator 
in a Banach space X , such that its resolvent set contains (—oo,0) and the resolvent satisfies: 

||A(A + A)" 1 1| <M, A > 0. 


Then, for rj £ (0,1), 

(4.10) A'# = Sm(7n?) I™ + Ay^ds. 

77 J o 


Note that there is another definition, through singular integrals, of the fractional powers of the (—A)~ a 
in L p (R n ) frequently used in the literature. See [28] Chapter 2.2] for the proof of equivalence of the two 
definitions for 1 < p < see also [TT] section 4.3]. 

Moment inequality. We recall here, the moment inequality estimate valid for fractional powers of 
non-negative operators. They are suitable to compare various fractional powers. Recalling [381 p. 98], 
we have the following version of the moment inequality with precise constant; for 0<a</3<7<l, 


II^Vll < 


(sin 7r( 7 /3 _ a a) )(7-q) 2 

7r(7-0)OS-a) 


(M- 




where tp £ D^A 1 ). 


4.3. Moment inequality extended. The task here is to extend the moment inequality to the form 
suitable to compare the powers 1 + a and 1 (here a > 0 near 0 + ). We need to use a more general than 
(14.101) expression (11.61) for the fractional powers, taken from [28] (3.4), p.59], which states that: 

(4.11) A 1+a (p = ^ Sm(7r . Q) r A“[A(A + A)~ 1 ] 2 cj)dX, 0 £ D(A 2 ), 

(1 - ot)i r J o 

r( 2 ) 

where the original term r(a)r( 2 -a) bas been transformed using the known properties of the T function; 

r(l + o) = oT(a), T(a)T(l — a) = -r-r —r, T(2) = 2. 

sm(7raj 

We are using the following bound, valid for positive operators , in the calculations below: 

(4.12) ||[Al(A + Al]- 1 ||. Y <M + 1, 


splitting the integral over (0 ,L) and (L,oo), we estimate the first integral, 

(4.13) 


2sm(7ra) f L A „^ A + = 2sm(7ra) f L X <*( X + A yi[ A ( X + A yi} A(/)dx 


(1 — a)i 


(1 — a) 7 


2 sin(7ro:) 

(1 - u)tt J Q 


j xa YTx(M + r,d\\\A < p\\ x . 


We thus have: 
(4.14) 


. 2 sin(7ra) f L , ll9 , „ ,, 2 sin(7ra) „ ,., , . L a ,, . ,.. 

1 K J I \ a [A(X + A^fcjxiXWx < M[M + 1) — ||^||x- 


(1 — a)i 


(1 — a)7 
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The second integral over (L, oo) is estimated next: 
(4.15) 


2 sin(7ra) f°° .. i, 2 . 2 ,,,,, 2sin(7ra) f°° / M N 2 II . 2 , „ ,, 

-'— 1 y L A + } A 2 <j)d\\\x < J J L A (—) IIAVlkdA 


(l-a)7r 


2 sin(7ra) M 2 


A 


-L“- 1 ||A 2 0|| a -. 


(1 — a) 7r 1 — a 

Minimizing with respect to L > 0 we get L m i n = , which leads to the final estimate: 


(4.16) 


||A 1+C V \\x < 2Sm{lT ^ M(M + \)Lkl I 
(1 — a)i t a 


x + 


2 sin(7ra) M 2 
(1 — a)ir 1 — a 


JwPV llx 


2sm(y) M i +a(M + i )Q H^ll^ll^^a ^ 2M||^|| 

(1 — a) 1 no. 


x , 


as a 


0+. 


5. Conclusion. 

The PROBLEM connected with the 3-D Navier-Stokes equation is not solved in the paper. But we 
are explaining its reason. The viscosity in the classical N-S equation is too weak to prevent better norms 
of its solutions from blowing-up in a finite time, especially in 3-D case. The 2-D equation is still critical 
with respect to the standard L 2 estimates. In order, for small initial data we enjoy the standard property 
of problems with quadratic nonlinearity , that near zero function the estimates of solutions obtained from 
the equation are proportional to the square of the norm, therefore there is a ball centered at zero the 
solutions originated in it will never leave that ball. For such solutions, the nonlinear term in N-S is 
subordinated to the main part operator. Concluding; to solve the problem with the N-S one need to find a 
better a priori estimate than the standard one in L 2 , coming from the ’symmetry ’ of that equation. Then, 
the nonlinear term would be subordinated to the main part operator preventing possible blow-up of better 
norms of the solutions, and the global in time existence-uniqueness theory will apply. In particular, if the 
estimate in L N (fl) (L N+ (f2)) is known, it will make the N-D Navier-Stokes equation critical (sub-critical) 
due to the following calculation 

(5.1) \\P(u-V)u\\ x _i < c|||M| 2 ||[ i 2 (Q)]N- < C , ||w||^ 7 2^_ (j2)]jv ||u||[ L Af(n)]JV < c"||ti|| x i ||u||[z,W(f2)]JV. 

It seems that the fractal regularizations; m in 2-D and (12.611 . s > y, in 3-D, having much better 
properties of solutions, offer an alternative way of description of the flow. But until very recently nobody 
was using fractal equations as models of physical phenomena. 
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